Wetting on nanorough surfaces 
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We present in this letter a free-energy approach to the dynamics of a fluid near a nano-structured 
surface. The model accounts both for the static phase equilibrium in the vicinity of the surface 
(wetting angles, Cassie-Wenzel transition) and the dynamical properties like liquid slippage at the 
boundary. This method bridges the gap between phenomenological phase-field approaches and more 
macroscopic lattice-Boltzmann models. 



Viscous dissipation is a major problem in micro or 
even nanofluidic systems Large pressure gradients 
are indeed necessary to produce a significant flow at small 
scale. This effect is well illustrated by the Poiseuille law 
relating the flux to the size of the duct: for a cylindrical 
duct of radius R, one can check that the flux decreases 
like R 4 when a given pressure gradient is applied. The 
Poiseuille law is however obtained by assuming no-slip 
boundary conditions (BC) for the liquid at the walls, 
and recent experiments have shown that this might not 
always be the case 0]. In the presence of partial slip BC, 
one can expect the situation to be less penalizing, spe- 
cially if the slip length associated to the boundary is of 
the order of the duct radius or even larger 0]. Super- 
hydrophobic surfaces have thus attracted a considerable 
interest recently as potential candidates for highly slip- 
ping surfaces: if the liquid is repelled by the surface, one 
can imagine that a thin layer of gas (air or vapor) at the 
surface could produce a lubricating effect. The situation 
is however more complex since natural or artificial super- 
hydrophobic surfaces are extremely rough: roughness is 
indeed a key ingredient in superhydrophobicity, since it 
favors the trapping of air or vapor bubbles at the bound- 
aries 0| • This effect is due to capillarity, and the question 
to understand the coupling between the three ingredi- 
ents, hydrodynamics, capillarity and surface roughness 
is a true challenge. From the experimental point of view, 
the situation is still very controversial, with measured 
slip lengths varying on several orders of magnitude 
It is thus crucial to have numerical models able to inves- 
tigate the static and the dynamic properties of a fluid at 
a structured boundary. 

Several approaches have been considered (see [3| for 
a recent review), ranging from molecular dynamics at 
the nanometric scales [j| to hydrodynamic models (in- 
cludingthe lattice-Boltzmann approach 0) at micronic 
scales 
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and a very recent phase-field 
model at intermediate scales [T3|. The method we shall 
consider in this letter is intermediate between phase- 
field phenomenological approaches, not accounting for 
the specificities of the gas phase (low viscosity and high 
compressibility) and the more macroscopic lattice Boltz- 
mann approaches. It accounts explicitly for a possible 
liquid-vapor coexistence in the vicinity of a structured 
surface, and also accounts for an hydrodynamic flow with 
liquid slippage at the surface. This model is able to re- 



produce both static properties such as wetting angles or 
Cassie-Wenzel states, but also to predict the effective slip 
lengths of structured interfaces by accounting explicitly 
for the low viscosity gas layer that forms at the bound- 
ary. Moreover, this approach gives access to the intrusion 
dynamics of the fluid inside the pores of the structured 
surface, which is an important phenomenon when the 
liquid is submitted to large pressures. 
The model is constructed as follows: 
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is the transport equation for the local density field p(r, t), 
where t is time and r denotes the position; v(r,t) is 
the local velocity field and Q[p] = F[p] — p J p(r)dr is 
the grand-potential functional associated to the thermo- 
dynamic equilibrium between the liquid and the vapor 
phase (F[p] is the free-energy), p is thus a physical order 
parameter here, and not simply a mathematical object 
to identify the phases as in usual phase-field approaches 
The convective transport of the density field is in- 
cluded in the left term, while the local thermodynamic 
equilibrium is contained in the right terms: the first writ- 
ing corresponds to the Allen-Cahn model interest- 
ing for the study of the phase diagram (non-conserved 
dynamics, fixed chemical potential p), and the second 
writting is the Cahn-Hilliard model [la |. corresponding 
to a more realistic conserved dynamics, that we shall use 
to investigate the non-equilibrium properties. C is re- 
lated to the molecular diffusion constant I? in a bulk 
phase by the relation: D — C'f"(p) where /" denotes 
the second derivative of the bulk free energy f(p) at the 
bulk density p. In this work we used D = 2300pm 2 /s 
in the liquid phase, to match the auto-diffusion con- 
stant in pure water (D — 2270pm 2 / s). The free energy 

functional is F[p] = J dr {f(p) + ^|Vp| 2 + pV waU (r)} , 

with f(p) = k B Tp {log (i^) - l} - ap 2 the free en- 
ergy density of the bulk van der Waals theory (a is the 
mean field attractive energy between two atoms, and 
b is the close packing inverse density) . The square 
gradient term accounts for the natural thickness of the 
liquid-vapor interfa ce and gives rise t o the surface ten- 
sion 7lv = w f£ h y / 2(f(p) — pp + P)dp [16| where P is 
the bulk pressure, pl and p\ are the coexisting liquid 
and vapor densities respectively, w is related to the in- 



2 



terfacial thickness. V wa n(r) is the interaction potential 
applied by the wall. Equation (UJ) has to be solved with 
the corresponding equation for the velocity field: 
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dt 
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where p m {r,t) is the local mass density (p m (r,t) = 
Mp(r,t), with M the molecular mass), a = rj(p)(\ / v + 
Vi>*) is the local viscous stress tensor (with i](p) the 
shear viscosity). Although more accurate prescriptions 
can be considered, we use in this work the simple ansatz 
7 l(p) = Pm^, where the kinematic viscosity v is assumed 
to be the same in the liquid and the vapor phase. The sec- 
ond term — pV ^ SF $f i ) is the thermodynamic force field 

applied by the density field p(r) to the flow. This term 
contains both the compressibility of the fluid, fixed by the 
van der Waals theory, and the capillary force when inter- 
faces are present. The last term of equation ([2]) accounts 
for the interaction with the wall. 

We shall now discuss the interactions between the fluid 
and the wall. The two fields p(r,t) and v(r,t) are de- 
fined in the whole parallelepipedic resolution box, and 
are thus defined and solved inside the walls. While the 
static properties of the wall are controlled by V^i^r) 
(expulsion of the fluid from the wall and wetting prop- 
erties), the boundary conditions for the velocity field are 
controlled by / wa ii. A prerequisite for a correct descrip- 
tion of the dynamical behavior of the liquid at super- 
hydrophobic surfaces is to model properly the intrinsic 
slip properties of the liquid on a flat solid surface. Our 
prescription for / wa n ensures partial slip BC in this case. 
These boundary conditions correspond to the continuity 
of the tangential stress criu at the fluid/solid interface: 
<7iu in the fluid equals — Av\\(z s ), the friction stress ap- 
plied by the wall, where z is the normal direction, || the 
tangential direction under consideration, A is the friction 
coefficient and vn (z s ) the slip velocity at the wall (z = z s ). 
The partial slip BC thus expresses: 
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where b = i]/X is the slip length, and r] is the shear vis- 
cosity of the fluid. If b, r\ and the slip velocity are known, 
the stress at the interface is —r]v\\(z s )/b. 

The bare interaction between an atom of the wall and 
an atom of the fluid is the Lennard- Jones potential in 
our model: Vtj(r) = 4e [{cr/r) 12 - (a/rf] . The global 
potential applied by the wall on a single atom of the fluid 
is V^uir) = f /? wa u(T-')^(k' - r\)dr', where p wan {r) is 
the density of atoms in the wall (assumed to be a con- 
stant inside the wall and zero elsewhere), and V(r) is 
the bare interaction potential [VljM] multiplied by the 
pair distribution function between an atom of the wall 
and an atom of the fluid. The main property of this 
pair distribution function is to vanish when r — > since 
two atoms cannot overlap, while VLj(r) — > +oo in this 
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Figure 1: Schematic views of the systems under consideration. 



limit. The product of these two functions thus goes to 
a well defined limit in r = 0. As a result, the poten- 
tial energy inside the wall is large, but finite. We shall 
denote by Vcut this energy and use the simple ansatz: 
F wa ii(r) = min(V r cut ,J p wa n{r')V LJ (\r'-r \)dr'). This pre- 
scription fixes the interaction to be of the Lennard- Jones 
type close to the boundary with the fluid, and to be V cut 
inside the wall. V cut fixes the value of the fluid density 
inside the wall, chosen to be negligible as compared to 
the density of the vapor phase (pa 3 < 1CP 3 ). In this 
limit, the results become independent of V cut - 

The force field / wa ii exerted by the wall on the fluid is 
chosen as a friction force: 



wall 



kep(r)v(r) / dr' < p wan (r')e 
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. This force is proportional to the interaction energy e 
between the solid and the fluid 0, [jjj ■ This prescription 
enabled us to reproduce correctly the evolution of the BC 
with the wetting properties. In particular, we obtain a 
transition from a no-slip BC in a wetting situation (9 ~ 
°) to a partial slip BC with b ~ 10 nm (depending on 
the value of k) in a non- wetting situation (9 ~ 120 °). 

The two equations ^ and ^ can be solved numeri- 
cally on a regular lattice (cubic unit cell in 3D or square 
in 2D) with periodic BC at the edges of the resolution 
box on the appropriate quantities, as we shall discuss. 
The wall is defined in the resolution box by the function 
Pwaii( r )- We consider a slab geometry corresponding to 
a fluid confined between two walls, the bottom wall is 
textured with a periodic structure (crenels in this study), 
while the top wall can either be structured in static stud- 
ies (FigJTJc), or planar in dynamic situations (Fig[Tj-d). 
The density p(r) is periodic in the slab (x-y) direction, 
but also in the z direction since it goes to a constant value 
(nearly vanishing) inside the walls. The velocity field is 
periodic in the x-y direction, but not in the z direction 
where a linear shear v s (r) is applied. In this case, the 
periodic BC is taken on u(r) = v(r) — v s (r) in any di- 
rections. By using a simple Euler scheme, equations ^ 
and (121) can be solved iteratively starting from an initial 
configuration for p(r) and v(r) . Taking advantage of the 
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periodicity of p(r) and u(r), we use an implicit method 
in the Fourier space to improve the numerical stability, 
but simpler schemes can be used as well. 

To illustrate the ability of the model to account for the 
small scale physics of the interfaces, we consider two situ- 
ations. The first study corresponds to a static case, where 
v s (r) — 0. The wetting properties can be probed ei- 
ther by measuring the contact angles of a drop placed on 
the surface, or more simply by measuring the three sur- 
face tensions of the liquid-solid (7sl)> vapor-solid (7sv) 
and liquid- vapor (tlv) interfaces. The model presented 
here gives a direct access to the free-energy or the grand- 
potential (an advantage compared to molecular simula- 
tions) ; the surface tensions can thus be measured directly 
by considering a fluid confined between two structured 
walls as depicted in FigQJc. Provided the distance be- 
tween the walls is large compared to the typical relax- 
ation length of the density profile, the walls can be con- 
sidered as independent. The relaxation length scale of 
the density profile is given by the interfacial thickness w, 
and we use a distance between the two walls of the order 
of 50u> , which ensures an accurate determination of the 
surface tensions. The liquid-solid [resp. vapor-solid] sur- 
face tension is measured by starting the simulation with 
a liquid [resp. vapor] phase between the walls. For planar 
interfaces this prescription is very efficient and the two 
surface tensions can be measured in the full range of vari- 
ation of the interaction potential Vlj , which is done by 
varying e. From these two quantities, and the knowledge 
of 7lv, the contact angle with the planar surface can be 
calculated through the quantity T = 7s ^~J SL . T corre- 
sponds to cos (9 when |T| < 1 and is an increasing func- 
tion of e (attractions favor wetting of the liquid phase) . 
We shall consider T rather than e in the discussions, since 
the contact angle is the experimentally relevant quantity. 
When the wall is not planar, and has a crenel shape as 
depicted in FigQJa, the situation is more complex since 
various states can be observed. We used in this case 
two different types of initial configurations: Wenzel states 
corresponding to a monophasic system, or Cassie states 
corresponding to a diphasic situation (see FigQ]-c). Of 
course, the equilibrium result should not depend on the 
initial state, this prescription gives us however a way to 
probe metastability in the system. The wetting prop- 
erties of the structured wall can be measured as well by 
determining the effective surface tensions 7gy and 7£y for 
the equilibrated configurations, the variation of the corre- 
sponding wetting parameter r eff is plotted in FigO This 
figure illustrates very well the Cassie- Wenzel transition 
between the imbibed and super- hydrophobic states Q. 
Moreover, several important features can be noted from 
these results: first we observe a strong metastability, with 
the presence of super hydrophobic metastable states, and 
next, we observe an asymmetry in the wetting properties 
when we go from the non-wetting region to the wetting 
one. The macroscopic model describing the transition 
between Cassie and Wenzel states is compared to the nu- 
merical results in Fig|2l This model is quite accurate in 
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Figure 2: Wetting properties of a crenelated surface (see text 
for details). 

the non-wetting region, which is not a surprise since we 
have the same geometry: the interfaces are planar at co- 
existence. In the wetting regime on the contrary, we ob- 
serve discrepancies with the macroscopic theory. The dif- 
ference comes from corner effects: the interaction poten- 
tial between the fluid and the wall has particular values 
at the corners that generate a corner energy. These cor- 
ner effects are mainly visible when a liquid phase fills the 
crenels (wetting situation), in gas phases the low value 
of the density leads to negligible contributions. If eu and 
cd are the excess free-energies in the liquid phase at the 
upper and the lower corners (see in Fig.l-c the points U 
and D), the Wenzel theory can be modified as follows: 

r eff = r(i + 2d/ l) - 2 £u + £p (4) 

where eu and cd are measured independently. This 
theory reproduces very well the effective contact angle 
(Figl21). Please note that eu and cd are not constant, 
they vary slightly with cos 9, which is visible in Fig|2l 

We now turn to the dynamical properties of the sur- 
faces. The dynamics of a fluid close to a textured wall 
can be accounted for using effective partial slip BC 
where both the slip length b and the hydrodynamic po- 
sition of the wall z s are unknown. For a smooth sur- 
face, it has been shown that the BC applies at about one 
molecular layer inside the fluid, for both no-slip and slip 
cases pjj. However for superhydrophobic surfaces, it is 
not clear where the BC should apply. In this work we 
therefore measured both b and z s , by probin g ou r sys- 
tem using both Couette and Poiseuille flows [12). The 
geometry of the system corresponds to FigQJd), where 
the top wall is planar and the bottom one is structured. 
This hybrid geometry allows for the determination of the 
effective boundaries for both the planar and the textured 
interfaces at the same time. The figure corresponds to 
a Couette flow, for which the top wall moves with a ve- 
locity U perpendicular to the crenels, while the bottom 
wall is fixed. The Poiseuille situation corresponds to fixed 
walls, but a pressure gradient is applied. We shall inves- 



4 



Both graphs: 

A A flat wall 

O-O SH surface 



% 

AAA&AAA A A A - A — A--- A 



5 

8 (nm) 



A-A-A-A -A-A A A A 

Periodicity L=30nm "q 
Width l=27nm 
Depth d=30nm 



5 
S (nm) 



Figure 3: Evolution of the slip length b with the recess 5 of 
the meniscus inside the crenel, for a flow perpendicular to the 
crenel. Inset: Evolution of the shear plane position z s as a 
function of 6. For the SH surface, z w is defined as the top of 
the crenels. 



tigate the variation of both b and z s as a function of 
the liquid pressure. Rather than the pressure itself, we 
can equivalently fix the average density in the slab by 
using a conserved dynamics, and measure the pressure 
afterwards. We shall use 8 , the recess length of the liq- 
uid in the crenel (taken from the top of the crenels, as 
depicted in FigQJd), as the control parameter. 8 can 
either be negative or positive. The slip length of the 



planar (bare) liquid-wall interface is taken to be around 
13 nm at coexistence, and is observed to vary very weakly 
with pressure (Fig[3|). On the contrary, the slip length of 
the structured wall varies in larger proportions, between 
35 nm for the upper limit down to for the lower one. 
Pressure has the effect to reduce slippage, and even to 
cancel it. The effective slip length crosses the bare one 
when 8 C ~ 2 nm, and is lower above. This result is largely 
independent of the crenel depth since the depth we con- 
sider in the present simulation (30 nm) is large compared 
to the crossover length S c . 



To conclude, the model presented here is very flexi- 
ble and can account for complex BC. Although we fo- 
cussed on crenelated walls, we can treat any type of ge- 
ometry, 3D shapes such as posts or random surfaces, or 
conic shapes to model surfaces coated with nanotubes 
[l8l |. The dynamics on crenelated surfaces revealed the 
strong sensitivity of the effective slip length to the recess 
of the meniscus 8 inside the crenels, cancelling slippage 
as soon as 8 ~ 2nm. This critical value increases with 
the periodicity L, although it remains a small fraction of 
it. However, increasing L is not necessarily a good idea 
since the corresponding intrusion pressure will decrease. 
Although nanostructured surfaces have lower slip lengths 
than microstructured ones in general, the slip properties 
resist better to pressure. There is thus a compromise to 
find between slippage and pressure resistance. 
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